Survival of parity effects in superconducting grains at finite temperature 
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We study the thermodynamics of a small, isolated superconducting grain using a recently developed quantum 
Monte Carlo method. This method allows us to simulate grains at any finite temperature and with any level 
spacing in an exact way. We focus on the pairing energy, pairing gap, condensation energy, heat capacity and 
spin susceptibility to describe the grain. We discuss the interplay between finite size (mesoscopic system), 
pairing correlations and temperature in full detail. 
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The bulk properties of a superconductor are well described 
by standard BCS theory. When the system size is reduced 
however, its mesoscopic behavior is strongly dictated by the 
finite electron number. For such small systems with a fixed 
number of particles, BCS theory is no longer applicable since 
the BCS order parameter is identically zero. Therefore it can- 
not determine the lower size limit for which the system ex- 
hibits superconducting properties. It was suggested by Ander- 
son 1 1 ] that superconductivity would disappear once the aver- 
age level spacing d (« 1 /V, V being the volume of the system) 
of the electron spectrum becomes larger than the bulk super- 
conducting gap A. Due to a series of experiments by Ralph, 
Black and Tinkham (RBT) 1 2] on the transport through a sin- 
gle superconducting nm-scale Al grain, a lot of authors shed 
new light on Anderson's suggestion. In their experiments, 
RBT found a spectroscopic gap larger than the average level 
spacing, which goes to zero when applying a suitable mag- 
netic field. The measurements also revealed a peculiar parity 
effect: grains with an even number of electrons have a larger 
gap in the spectrum than grains with an odd electron number. 
These observations were regarded as signs of 'superconduc- 
tivity', in the sense that there is a pair-correlated ground state. 
Properties indicative of strong pairing correlations were only 
found in grains with d < A. So Anderson's answer turned out 
to be uncomplete, since it does not differentiate between odd 
and even numbers of electrons. A large number of theoreti- 
cal studies tried to characterize the ground state correlations 
and superconductivity of such small systems in a qualitative 
way and tried to predict the critical level spacing at which the 
superconductivity breaks down. An extended review can be 
found in Ref. |3[. In this report we study the competition 
between pairing, finite size and finite temperature in an exact 
way, with all quantum correlations taken into account. 

To model small superconducting grains, one uses the re- 
duced BCS Hamiltonian 0|: 

a a 
H= £ (Ej-on B h)cl c c La -Xd £ B)B f , (1) 

o=±,;=l j,j'=l 

where = ct + cj_. The operator cj a creates an electron 
in the single-particle state \j,o). The quantum number j la- 



bels the Q. single particle levels with energies £,-, and o labels 
time reversed states. Since the pairing interaction only scat- 
ters time-reversed pairs of electrons within an energy cod of 
the Fermi level £f, electrons outside the cutoff are not taken 
into account. The parameter A, is the dimensionless BCS cou- 
pling constant and is related to A and cod via the bulk gap 
equation sinh(l /X) = <Md/A |5]. We take X — 0.224, close to 
that of Al 1 6]. The Zeeman term couples an external magnetic 
field h to the electrons and /jg is the Bohr magneton. Through- 
out the paper, we will consider a half-filled band with fixed 
width 2(0o and Q, = 2<x>o jd doubly degenerate and uniformly 
spaced levels with energies Zj — jd. We will only discuss the 
case without magnetic field h. 

To study the cross-over from the bulk to the few electron 
limit, a number of authors originally used a parity-pro j ected 
grand canonical (g.c.) BCS approach HHUBUElIj- The 
parity effect can be explained with this variational technique. 
However, an artificial sharp transition to the normal state ap- 
pears at some critical level spacing and temperature, which 
is impossible for a finite system. Since the electron number 
fluctuations are strongly suppressed by charging effects in the 
experiments of RBT, it is clear that a canonical formalism is 
needed to describe the grains properly. A number of canon- 
ical techniques were used to tackle this problem. Unfortu- 
nately, exact diagonalization techniques (e.g. Lanczos lol l 
can only handle systems with a very small number of elec- 
trons. In order to go to larger model spaces, particle num- 
ber projection was combined with the static path approxima- 
tion (SPA) plus random-phase approximation (RPA) treatment 
HI 3l 1 1411 and with variational wavefunctions |6]. Dukelsky 
and Sierra developed a particle-hole version of the density- 
matrix renormalization-group (DMRG) method to study the 
crossover I15lllal . All these canonical techniques reveal the 
parity effect at low enough temperatures, and make clear that 
the abrupt cross-over is just an artefact of the g.c. approach. 
It turned out that small grains with d < A are indeed char- 
acterized by strong superconducting pairing correlations. As 
the grain size decreases, quantum fluctuations of the order pa- 
rameter start to play a crucial role. These fluctuations make 
the cross-over completely smooth without any sign of critical 
level spacing. Only when the grain is not too small (d <C A) 
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the fluctuations in the order parameter can be neglected, mak- 
ing the mean field description of superconductivity appropri- 
ate. In the canonical picture, pairing correlations still exist 
at arbitrary large values of d/A, though in the form of weak 
fluctuations. Qualitative differences between the pairing cor- 
relations in the bulk and the few-electron regime make it still 
possible to speak of the superconducting regime (d <C A) and 
the fluctuation-dominated regime (d > A) yj. 

It was only after the appearance of most of these works that 
one became aware of the fact that the reduced BCS model has 
an exact solution, worked out decades ago by Richardson in 
the context of nuclear physics 1 1711 . In Ref. Ill 811 . Sierra et al. 
compare the previously mentioned treatments with the exact 
solution. Using this exact solution to study the finite tem- 
perature behavior for a large number of many-particle states 
is difficult due to the exponential scaling of the number of 
eigenstates that need to be considered. Gladilin et al. devel- 
oped an approximation based on the Richardson solution to 
get finite temperature information 1 19]. In Ref. |3J it was al- 
ready suggested by von Delft and Ralph that quantum Monte 
Carlo (QMC) techniques could be helpful to investigate the 
BCS pairing model at finite temperature. Recently we devel- 
oped a new quantum Monte Carlo method 12(1 \2 11 that is able 
to simulate the BCS model for any fixed number of particles 
without sign problem. The method allows calculating ther- 
modynamic properties in an exact way, up to a controllable 
statistical error. Simulations can be performed at any finite 
temperature and any level spacing d/A for large system sizes. 
Because our method allows a projection on specific symme- 
tries like the total spin projection, we can calculate the sus- 
ceptibility and magnetization. 
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FIG. 1 : The temperature dependence of the pairing energy per elec- 
tron for grains with an even number N = £2 of electrons. Simulations 
were performed for different grain sizes £2. The energy and temper- 
ature scale is set by the level spacing d. 

We performed simulations of grains with different sizes (£2 
equal to 10, 40, 80, 160, 320 and 400). These half-filled 
model spaces lead respectively to ratios d/A of 8.68, 2.17, 
1.09, 0.54, 0.27 and 0.22. Figs. [j]and|2]show the thermal av- 
erages of the pairing energy Hp — — 'kdYffy=\ "j^f P er parti- 
cle as a function of temperature for even and odd grains. The 
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FIG. 2: The pairing energy per electron as a function of temperature 
for an odd grain with different sizes. Simulations were performed for 
grains with N = £2 + 1 electrons. 



energy scale is set by the level spacing d. By comparing both 
figures, one notices that at low enough temperatures (typically 
T <d) the even electron system has more pairing energy than 
the odd system. This is due to the single unpaired electron, 
which blocks the Fermi level in the odd case. Around T ~d a 
small dip appears in the odd pairing energy. Qualitatively, this 
can be explained as follows: due to the thermal energy, the 
single unpaired electron is moved one level upward, making 
the Fermi level available to pair scattering. This is reflected in 
a slight decrease of the pairing energy in Fig. [2] To measure 
the real correlation energy due to pairing in the system, the 
'canonical' pairing gap 

AL = (^O 2 t « B ™ B »> - <4»»>JW>), ( 2 ) 
m,n=l 

was introduced in Eq. (92) of Ref. 0. The second term 
subtracts the thermal average of the pairing interaction for the 
non-interacting system. When going to the thermodynamic 
limit, Acan becomes equivalent to the BCS bulk gap A jj]. 
Figure [3] shows the even and odd canonical gap for differ- 
ent system sizes. It follows very clearly that the temperature 
scale at which the parity effect appears is set by the level spac- 
ing d, and this for all grain sizes. The crossover temperature 
is given by T cr = AlnAf e ff, with A^ff the effective number of 
states available for excitation (N e ff = V8%TA/d in the limit 
d <C A) I22I1 . This is in qualitative agreement with Figure [3] 
where the crossover temperature decreases as the grain size 
is reduced. One should of course keep in mind that the tem- 
perature is shown in units of the level spacing which is con- 
siderably smaller for the largest grains. Figure shows that 
pairing correlations persist even for ultrasmall grains and that 
a reduction of the grain size leads to a suppression of these 
correlations. 

The condensation energy E con & = (\\t\H\\\t) — (FS\H\FS) is 
the energy difference of the state \\\t), where all quantum cor- 
relations are included, and the non-correlated Fermi sea \FS) . 
Figure |4] shows the thermal average of the condensation en- 
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FIG. 3: The canonical pairing gap as a function of temperature. For 
each number of levels £1 the gap is calculated for an even (N = £1) 
and an odd (N = £1 + 1) number of electrons. Only at low enough 
temperature one can distinguish between the gap of the even grain 
(upper curve) and the odd grain (lower curve) of the same size £1. 




FIG. 4: The condensation energy per particle as a function of T/d 
for system sizes £1 = 10,80 and 400. Even (odd) grain data points 
are connected by a solid (dashed) line. 



ergy per particle for a number of even and odd grains as a 
function of temperature. These energy differences were ob- 
tained by calculating the thermal averages of the Hamiltonian 
over correlated states \\\t) and over the Fermi states \FS) sepa- 
rately. Below temperatures of the order d, the even grains have 
a larger condensation energy (in absolute value). Both even 
and odd grains have a minimal condensation energy around 
T w d. In agreement with Ref. JjJ], our calculations give a 
quasi intensive condensation energy for the smallest grains 
(d ^> A), while the condensation energy of grains with d <C A 
increases (in absolute value) inversely proportional to d. 

Figure|5]shows the heat capacity as a function of tempera- 
ture for sizes £1 = 10, 80 and 400. Around the crossover tem- 
perature where the parity effect becomes visible (see Figures 
[5]and|4j, a slight parity effect also appears in the heat ca- 
pacity. Here the even heat capacity exceeds the odd one. At 
higher temperatures the odd and even results become indistin- 



guishable again. For the £1 = 10 grain size, the finite model 
space makes the Shottky peak visible when the temperature 
becomes of the order of the level spacing. 




FIG. 5: The heat capacity c = -jS)- as a function of T/d for system 
sizes £1 = 10, 80 and 400. Even (odd) grain data points are connected 
by a solid (dashed) line. Around temperatures T ~ 0.5rf for £1 = 
10, 80 and T ~ d for £1 = 400 the heat capacity of the even grain (with 
N = £1 electrons) exceeds the odd (N = £1+ 1) specific capacity. 



The spin susceptibility of a grain is defined by 



X(T) = 



dh 2 



h=0 
2\ 



(M) 



(3) 



Here J = — TlnZ is the free energy of the grain, with 
Z the canonical partition function. The susceptibility is 
proportional to the fluctuation of the 'magnetization' M — 
~f J BY l a,nG c n a c n,a at finite temperature T. The spin suscepti- 
bility of a single isolated grain has been studied by Di Lorenzo 
et al. l23ll . They found that the pairing correlations affect the 
temperature dependence of the spin susceptibility. In partic- 
ular, if the number of electrons in the grain is odd, the spin 
susceptibility shows a re-entrant behavior as a function of T 
for any value of the ratio d/A. They show that this behav- 
ior persists even in the case of ultrasmall grains, where the 
level spacing is much larger than the BCS gap. Since this 
re-entrance behavior is absent in normal metallic grains, they 
suggested that this quantity can be measured and used as a 
unique signature of pairing correlations in small and ultra- 
small grains. The susceptibility was calculated by combining 
an analytic analysis in the limiting cases A^S> d and A <C d 
with a static path approximation for intermediate values. By 
means of exact canonical methods based on Richardson's so- 
lution, they also got exact results at low temperatures. With 
the aid of our QMC method, we are now able to solve the 
problem exactly for the whole temperature range. Figures [6] 
and[7]show the temperature dependence of the spin suscepti- 
bility for a number of even and odd grains, respectively. The 
susceptibility is normalized to its bulk high temperature value 
%p = 2/jg/d. Our results are completely in line with those of 
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Di Lorenzo et al. l23ll . At low temperatures the even suscep- 
tibility remains exponentially small, while for an odd grain 
the unpaired spin gives rise to an extra paramagnetic contri- 
bution to the spin susceptibility (% ~ 1*b/T). The minima in 
the odd spin susceptibilities coincide with a small increase of 
the pairing correlations (see Figs. |2]and[3}, with a mimimal 
condensation energy (see Fig. |4} and with a parity effect in 
the heat capacity (see Fig. [5}. For the smallest odd grain no 
re-entrant behavior is visible in Fig. This is an effect of the 
finite model space. If the BCS coupling constant is increased 
a little, a re-entrance effect appears also in this case. 
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FIG. 6: The spin susceptibility normalized to its bulk high tempera- 
ture value %P as a function of T/d for a number of even grains. Each 
grain contains N = CI electrons, with £2 the model space size. 
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FIG. 7: The spin susceptibility as a function of T/d for a number 
of odd grains (containing N = £2 + 1 electrons). The grain size is 
determined by the model space size £2. 

In conclusion, we solved the BCS pairing problem at finite 
temperature exactly via quantum Monte Carlo simulation. We 
studied odd and even grains with a large number of electrons 
and arbitrary level spacings. Our exact results confirm predic- 
tions of previous approximate calculations, showing that the 
physics of ultrasmall superconducting grains is well described 
by a pairing model with exact particle number projection and 
that parity effects are visible in thermodynamic properties. 



The number of unpaired electrons in a grain can be in- 
creased by an external magnetic field. Frauendorf et al. 
showed that at zero temperature a magnetic field attenuates 
the pairing, but for a mesoscopic system in a strong mag- 
netic field the pairing correlations may come back after heat- 
ing l24ll . Such a re-entrance of pairing correlations has also 
been discussed by Balian et al. Ill ill . Work on this problem of 
how an external field can influence the thermodynamic prop- 
erties of a single superconducting grain is in progress. 
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